Lattice artefacts and the running of the coupling constant. 
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We study the running of the Liischer-Weisz- Wolff (LWW) coupling constant in the two dimensional 
0(3) nonlinear o model. To investigate the continuum limit we refine the lattice spacing from the 
-L value used by LWW up to -^r-. We find much larger lattice artefacts than those estimated by 
LWW and that most likely the coupling constant runs slower than predicted by perturbation theory. 
A precise determination of the running in the continuum limit would require a controlled ansatz of 
extrapolation, which, we argue, is not presently available. 
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The hallmark of QCD is its alleged asymptotic free- 
dom (AF), that property which expresses the fact that 
at shorter distances the interactions between quarks and 
gluons are becoming weaker. This fact however has been 
established only in perturbation theory (PT), an approx- 
imation scheme without a mathematical basis and which, 
moreover, has been shown to be plagued by ambiguities 
(expectation values of variables of compact support de- 
pend upon the boundary conditions (b.c.) used to reach 
the thermodynamic limit) ]l|. 

It is therefore most important to establish whether AF 
is really a property of QCD in a nonperturbative frame- 
work. The first step in this direction was taken in 1991 
by Luscher, Weisz and Wolff (LWW) 0, who proposed 
a method to investigate the presnce of AF in the two 
dimensional (2D) O(N) nonlinear a models, which, per- 
turbatively, are also AF for N > 3. As a coupling con- 
stant they proposed the following renormalization group 
invariant: 



g 2 (P,L) = 



2L 



(A-1K(L) 



(1) 



Here £(£) is the correlation length of the O(N) model 
in an infinite strip of width L with periodic b.c. It is 
defined by the following double limit: consider a finite 
strip of size L x T with periodic b.c. in the direction of 
size L and arbitrary b.c. in the direction of size T. Then 



f(£) = -lim x ^ 00 lim T ^ c 



ln(< s(0) • s(x) >) 



For the standard O(N) action 



H, 



hi 



-s(i) ■ s(j) 



(2) 



(3) 



at inverse coupling constant j3, if in cqn.(l) one ex- 
pressed L in units of the thermodynamic correlation 
length £oo(/3)j then in the limit L — ► oo and j3 — > j3 cr t, 
holding z = L/£oo fixed, one would obtain a unique func- 
tion g 2 (z) describing the running of g 2 with the physical 
distance z. 



As we pointed out in 1992 ||, the interpretation of g 2 
as a coupling constant is somewhat misleading since it 
does not measure the strength of the interaction. Indeed 
in the the massive continuum limit of a free field theory 
g 2 (z) is a nontrivial function, running linearly with z. 
LWW argued in favour of their choice by pointing out 
that in PT, to lowest order in the bare (lattice) coupling 
constant g 2 (L) ~ A. This argument is also problematic 
since even if (3 crt — oo, to construct the continuum limit 
one would have to let L — > oo and thus reach a regime 
where PT in the bare coupling would clearly not be ap- 
plicable. 

Nevertheless g 2 (L, (3) is a renormalization group invari- 
ant and if one would discover that for some j3 < oo, g 2 (L) 
became independent of L for large L, that would mean 
that at that (3 the model is critical, which, as we will 
explain below, would rule out the existence of AF in the 
massive continuum limit. 

In 1991, for the 0(3) model, LWW claimed to be 
able to establish the continuum running of g 1 '(L / '£t»(/3)) 
up to physical distances as small as 0.03304(4) and to 
verify that it approached the perturbative (AF) predic- 
tion. To achieve this they employed a finite size scal- 
ing (FSS) technique: one measures £(L) at some j3 and 
L, then leaving (3 unchanged, one doubles L and mea- 
sures £(2X). One thus obtains a scaling curve giving 
g 2 (2L)/g 2 (L) versus g 2 (L). This step scaling curve, as 
LWW called it, allows one to connect small physical dis- 
tances to large ones (£/£oo (P) > 7) , where the continuum 
limit of g 2 (L/^oo(f3)) could be reached. 

Similar FSS approaches were used in 1993 by Kim Gj 
and in 1994 by Caracciolo et al fa] to predict the value 
of £(/?) up to £ ~ 10 5 , even though the largest lattices 
involved did not exceed L = 512. However whereas these 
authors produced their scaling curves simply by observ- 
ing that the Monte Carlo (MC) data coming from differ- 
ent values of /3 seemed to fall on the same curve, which 
they took as their step scaling curve, the LWW paper 
claimed to have really controlled lattice artefacts (the 
aproach to the continuum). More precisely the problem 
is this: of course if one knew the continuum value of the 



step scaling curve one could connect small physical dis- 
tances to large ones. But in an MC investigation, by ne- 
cessity, one can only gather data at finite cutoff (l/£oo(/3) 
or alternatively 1/L). However the continuum limit re- 
quires letting £00 (/3) — ■* 00 and L — > 00. Therefore one 
must in principle worry about extrapolating the results 
obtained for the step scaling function at finite cutoff to 
the continuum limit. 

This feat, which the other above quoted authors did 
not even attempt, was achieved by LWW by assuming 
a Symanzik type of approach to the continuum limit. 
Namely at fixed g 2 they assumed that the step scaling 
function approaches its continuum limit value as 1/L 2 
(strictly speaking, inspired by PT, the Symanzik fit in- 
volves log(L) J L 2 , however for 6 < L < 16 the log can be 
approximated by a constant). They backed this assump- 
tion with their MC data. However the values of L they 
used ranged only from 6 < L < 16. The main point of 
our paper is to show numerical evidence that if one goes 
to much larger L values, the Symanzik fit does not work 
and an entirely new picture for the running of g 2 with 
the physical distance emerges. 

Before showing our data, we would like to elaborate 
on a subtlety having to do with the numerical determi- 
nation of the correlation length £,(L). Namely while the 
definition in eq.(l) is mathematically well defined, one 
must adopt a computational procedure for implement- 
ing it. LWW used an L x T lattice with free b.c. in 
the T-direction, took T = 5L and assumed a pure ex- 
ponential decay for L < x < T — L (see Appendix B 
of ref . Q ) . We used instead the following numerical esti- 
mate of the correlation length £: let P — (p, 0), p 
n = 0,l,2,...,T-l, T = 10T. Then 



2n7T 

T 



i 



1 



2 sin(7r/T) 



y/G(0)/G(l) - 1 



(4) 



where 



G(P) = -^(\s(P)\ 2 ); s(P)="£e iP *s(x) (5) 



LT 



It is not clear whether the LWW prescription or the one 
adopted by us provides an estimate closer to the true ex- 
ponential correlation length. For finite T the LWW pro- 
cedure, employing free b.c, is likely to produce a value 
smaller than the true £(T). In the procedure adopted 
by us, there are two effects of opposite sign, which could 
bring the result closer to the true exponential correlation 
length: 

1. The periodic b.c. in the T-direction increases the or- 
der in the system compared to an infinite strip. 

2. Since our definition is sensitive to the multiparticle 
states, for T — > 00 it would produce a value smaller than 
the true £(£). This effect was studied by Campostrini et 
al B using the high temperature expansion and found to 
be less than 0.2%. 

Since we want to compare our data with those of LWW, 
we show in Fig.pl the resulting step scaling functions com- 
puted with the two procedures at the same (3 and L. As 



can be seen, within numerical accuracy, at these values 
of g 2 the two procedures produce similar results. 
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FIG. 1: Step scaling function g (2L)/g (L) versus g (L) for 
L — 20, 40, 80 computed with our method of estimating £(L) 
(solid lines) and with the LWW method (dotted line). 



1.2 



1. IE 



1.16 



531.14 



1.12 



1.1 



PT 



LWW continuum 
MC L=30 



0.8 



0.9 



FIG. 2: Step scaling function g 2 (2L)/g 2 (L) versus g 2 (L) for 
L — 20. Our MC data are connected by a spline. 

Our results for the step scaling function are shown in 
Fig.|, Fig.| and Fig.| for L = 20, 40 and 80 respectively. 
The values were obtained by measuring the correlation 
length at given (3 first on a lattice of size LxlOL, then on 
2i x 20L. The values of (3 used (1.815, 1.94, 2.05, 2.16, 
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FIG. 3: Step scaling function g (2L)/g (L) versus g (L) for 
L = 40. Our MC data arew connected by a spline. 
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FIG. 4: Step scaling function g 2 (2L)/g 2 (L) versus g 2 (L) for 
L = 80 and 160. Our MC data for L = 80 are connected by 



2.27, 2.38 and 2.49) were chosen so that g 2 took values 
in the range covered by LWW. The figures contain also 
two benchmarks taken from LWW Eqj: 

- The 3-loop PT curve. 

- The LWW estimated continuum values. 

Fig|| contains also our value for the step scaling func- 
tion obtained by doubling L from 160 to 320 at g 2 — 
1.05397(81). The error bars were estimated as follows: 
for each value of [3 and L we started from a randomly 



chosen configuration and ran the improved cluster algo- 
rithm using 100,000 clusters for thcrmalization and 
1,000,000 clusters for measurements. We then repeated 
this procedure a minimum of 157 times, except for the 
value at L — 320 which contains only 74 runs. We com- 
puted the average value over these samples of G(0) and 
G(l) and from them £(£). Since £(£) is a nonlinear func- 
tion of G(0) and G(l), we estimated the error for £,(L) 
and g 2 (L) from our sample of independent values by the 
jack-knife method. The data for g 2 (L) that were used in 
the figures are given in Tab.l. 

Our data do not agree with the LWW prediction g] 
for the step scaling function. The latter were obtained 
by extrapolating the MC values obtained by taking L 
between L = 6 and L = 16 via a second order polynomial 
in jp. We did not display the LWW values for 6 < L < 
16 because of lack of space. However when combined with 
our values, they reveal a rather complicated approach to 
the continuum. This is not an unexpected fact. Indeed if 
L is sufficiently small at given f3, the (asymmetric) lattice 
is strongly ordered in the transverse (shorter) direction. 
On the other hand, as we emphasized several years ago 
H , since in the continuum limit one must let L — > oo and 
the Mermin- Wagner theorem guarantees the restoration 
of the 0(N) symmetry in that limit for any finite (3, 
clearly for L sufficiently large this 'perturbative' regime, 
of spins highly ordered in the transverse direction, cannot 
persist. 

In fact we showed in [pj that even if (3 crt = oo, as 
predicted by PT, the spins would cease to be highly or- 
dered in the transverse direction for L sufficiently large. 
Indeed bare PT itself provides a clue as to the distance 
over which the spins remain well ordered since to lowest 
order one has 



MO) ■ *(*)) 



N - 1 



D(x) 



(6) 



and to a good approximation D(x) — \ + ^- ln(x). Thus 
bare PT suggests that spins are well ordered over dis- 
tances 0(exp([27r/3/(iV - 1)]). On the other hand the 
AF formula predicts £ = 0(exp[27r/3/(iV - 2)]). Thus at 
fixed physical distance (rh\), in taking the continuum 
limit one would surely leave the regime in which PT in 
the bare coupling is applicable. 

Returning now to the pattern of lattice artefacts, ini- 
tially, if (3 is large enough, at small enough L, they should 
follow the Symanzik pattern used by LWW because the 
system is essentially in a PT regime. This regime has 
nothing to do with the true approach to the continuum, 
which occurs only when L is sufficiently large (for given 
(3) so that the 0{N) symmetry becomes approximately 
true. How the true continuum limit is approached can 
be model dependent. For instance in the Ising model 
there are good reasons to expect a 1/L 2 leading behavior 
||. On the other hand in the 0(2) model there are both 
theoretical |Hj] and numerical reasons |JT[ to expect a 
l/ln(L) approach. 



It is reasonable to expect that the 0(3) model, enjoy- 
ing a continuous symmetry, behaves as the 0(2), not as 



the Ising model. We have attempted a 



ML) 



fit to our 



data but could not obtain a reliable continuum value. 
(We encountered a similar situation in investigating the 
renormalized coupling gn in the 0(3) model Q. The 
case of the 0(2) model is simpler to handle because the 
theory predicts both the continuum value and the leading 
correction |1 if . ) Even though we do not have a predic- 
tion for the continuum step scaling function, our results 
do not corroborate the original prediction of LWW and 
suggest that most likely the nonperturbative running of 
g 2 is slower than predicted by PT. This situation is con- 
sistent with, though in no way proving, the existence of 
a transition to a massless phase at finite f3 C rt, as argued 
by us recently 113]. In that paper we also proved rigor- 
ously that for the standard action, the massive contin- 
uum limit cannot be AF if j3 C rt < oo. The result follows 
from a Ward identity and the reflection positivity of the 
standard action. 

Finally, regarding the running of a s (Q) in QCD^, all 
we can say is that the Symanzik type of fit for the ap- 



J 



proach to the continuum has no justification there ei- 
ther. Indeed, that fit is inspired by PT. If in fact lat- 
tice QCD4 does undergo a deconfining zero temperature 
transition at nonzero (bare) coupling, so that the run- 
ning of a s (Q) does not follow PT, there is no reason to 
expect the lattice artefacts to follow the Symanzik ansatz. 
Therefore it would be very useful if the lattice community 
employed its resources to establish first the true cutoff ef- 
fects and the true running of a s in the pure Yang-Mills 
theory by going to larger L, before attempting to handle 
dynamical fermions; the latter unavoidably can only be 
done on minuscule lattices, and using the Symanzik fit to 
extrapolate to the continuum can be misleading, as we 
have found. As we said many years ago 1 14 , we expect 
that in the four dimensional Yang-Mills theory as well as 
in QCD4 there exists a nontrivial fixed point and that 
a s (Q) runs slower than predicted by PT, with the effect 
becoming pronounced by 1 TeV or less. 
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Tab.l: Monte Carlo data for g 2 (/? , L) 






1.815 


1.94 


2.05 


2.16 


2.27 


2.38 


L = 20 


1.06185(78) 


.91115(67) 


.81616(70) 


.74170(56) 






L = A0 


1.27230(85) 


1.04334(76) 


.91530(68) 


.81901(63) 


.74419(57) 




£ = 80 




1.24174(85) 


1.04778(72) 


.91907(70) 


.82404(60) 


.74811(53) 


L = 160 






1.24497(100) 


1.05398(81) 


.92457(69) 


.82564(58) 


L = 320 








1.2530(15) 
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